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,S^ • Abstract 

II I We have studied numerically the so-called Bondi-Hoyle (wind) accretion 

O ' on to a rotating (Kerr) black hole in general relativity. We have used the 

Kerr-Schild form of the Kerr metric, free of coordinate singularities at the 



C/3 ■ black hole horizon. The 'test-fluid' approximation has been adopted, assum- 



ing no dynamical evolution of the gravitational field. We have used a recent 
formulation of the general relativistic hydrodynamic equations which casts 
them into a first-order hyperbolic system of conservation laws. Our studies 
have been performed using a Godunov-type scheme based on Marquina's 
flux-formula. 

We find that regardless of the value of the black hole spin the final ac- 
cretion pattern is always stable, leading to constant accretion rates of mass 
and momentum. The flow is characterized by a strong tail shock which is 
increasingly wrapped around the central black hole as the hole angular mo- 
mentum increases. The rotation induced asymmetry in the pressure field 
implies that besides the well known drag, the black hole will experience also 
a lift normal to the flow direction. 

1 Introduction 

The term "wind" or hydrodynamic accretion refers to the capture of matter by a 
moving object under the effect of the underlying gravitational field. The canonical 
astrophysical scenario in which matter is accreted in such a non-spherical way 
was suggested originally by Bondi and Hoyle [pi, who studied, using Newtonian 



gravity, the accretion on to a gravitating point mass moving with constant velocity 
through a non-relativistic gas of uniform density. Such process apphes to describe 
mass transfer and accretion in compact X-ray binaries, in particular in the case in 
which the donor (giant) star lies inside its Roche lobe and loses mass via a stellar 
wind. This wind impacts on the orbiting compact star forming a bow-shaped 
shock front around it. 

The problem was first numerically investigated in the early 70's. Since then, 
contributions of a large number of authors using highly developed Godunov-type 
methods extended the simplified analytic models (see, e.g., |lj, |[ and references 
there in) . These Newtonian investigations helped develop a thorough understand- 
ing of the hydrodynamic accretion scenario, in its fully three-dimensional charac- 
ter, revealing the formation of accretion disks and the appearance of non-trivial 
phenomena such as shock waves or flip-flop instabilities. 

We have recently considered hydrodynamic accretion on to a moving black hole 
using relativistic gravity and the "test fluid" approximation [fl ^, [jj H ■ We present 
here a brief summary of the methodology and results of such simulations. We 
integrate the general relativistic hydrodynamic equations in the fixed background 
of the Kerr spacetime (including its non-rotating Schwarzschild limit) and neglect 
the self-gravity of the fluid as well as non-adiabatic processes such as viscosity or 
radiative transfer. In the black hole case the matter flows ultimately across the 
event horizon and becomes causally disconnected of distant observers . Near that 
region the problem is intrinsically relativistic and the gravitational accelerations 
significantly deviate from the Newtonian values. 

2 Equations 

The general relativistic hydrodynamic equations can be cast as a first-order flux- 
conservative system describing the conservation of mass, momentum and energy. 
Formulations of this sort are given, e.g. in |^, ^. In this work we follow the 
approach laid out in for a perfect fluid stress-energy tensor T^^'^ . The system 
of equation then reads: 



1 / d^u ^ d^^r 



(1) 



(x" = t; X* spatial coordinates, i — 1,2,3) where u = u(w) arc the evolved 
quantities, u — {D, Sj,t) and P are the fluxes 

r = {^D (.^ - ^) ,5, (.^ - ^) +p.;,. (.^ - ^) +p.^) , (2) 

v* being the 3-velocity and p the pressure. The corresponding sources s are given 

by 



We note the presence of geometric terms in the fluxes and sources which appear as 
the local conservation laws of the density current and stress-energy are expressed 
in terms of partial derivatives. These terms are the lapse function a, the shift 
vector /3* and the connection coefficients F^ of the 3+1 spacetime metric 

ds^ = g^.^dxt'dx'' = -(a^ _ p^p^)dt^ + 2P,dx'dt + j.jdx'dx^ (4) 

Additionally g = det{gfj,^) is such that V^S = Q^V^ and 7 = det(7ij). 

The vector w, representing the primitive variables, is given by w = {p,Vi,e) 
where p is the density and e the specific internal energy. The evolved quantities 
are defined in terms of the primitive variables as 13 = pW, Sj — phW'^Vj and 
r — phW^—p—D, H^ being the Lorentz factor H/^ = (1— -y^)"^/^, with u^ — ^ijV^v^ , 
and h the specific enthalpy, h = 1 + e + p/p. A perfect fluid equation of state 
p = (T — l)pe, r being the constant adiabatic index, closes the system. 

In our computations we specialize the above expressions to the Kerr line ele- 
ment which describes the exterior geometry of a rotating black hole. We use the 
Kerr-Schild form of the Kerr metric, which is free of coordinate singularities at 
the black hole horizon. Computations using the more standard Boyer-Lindquist 
(singular) form of the metric are presented in |p| . Pertinent technical details con- 
cerning the specific form of these metrics are given in pfl . 

3 Numerical scheme 

Our hydrodynamical code performs the numerical integration of system (|^) using 
a Godunov-type method. The time update from i" to t""*"^ proceeds according to 
the following algorithm in conservation form: 



u"+i = u 



At 



■^,J «J /S^^k 



(fi+l/2j - fj-l/2j) + At Si J , (5) 



improved with the use of (second-order) conservative Runge-Kutta sub-steps to 
gain accuracy in time |13| . The numerical fluxes are computed by means of Mar- 
quina's flux- formula [^ . After the update of the conserved quantities the primitive 
variables are computed via a root-finding procedure. 

The flux-formula makes use of the complete characteristic information of sys- 
tem (|^), eigenvalues (characteristic speeds) and right and left eigenvectors. Generic 
expressions are collected in [g|. 

The state variables, u, must be computed (reconstructed) at the left and right 
sides of a given interface, out of the cell-centered quantities, prior to compute the 
numerical fluxes. In relativistic hydrodynamics one has the freedom to reconstruct 
either w (primitive variables) or u (evolved variables). For efRciency and accuracy 
considerations we reconstruct the first set, from which the remaining variables are 
obtained algebraically. The code uses slope-limiter methods to construct second- 
order TVD schemes by means of monotonic piecewise linear reconstructions of the 
cell-centered quantities. We use the standard minmod slope which provides the 
desired second-order accuracy for smooth solutions, while still satisfying the TVD 
property. 
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Figure 1: Schematic representation of stationary supersonic wind accretion. The shock 
may be detached as in the figure (bow shock) or attached to the rear part (tail shock), 
depending on the flow asymptotic conditions (vao and Ca^) and thermodynamics of the 
gas (poo and V). 

4 Results 



The classical solution for an asymptotically uniform wind of presureless gas past 
a compact source (modeled analytically by a point mass) was obtained by |^. In 
this solution the material is focussed at the rear part of the object as a result 
of the gravitational pull. For a pressureless gas, the density at this symmetry 
line could reach an infinite value and matter would flow on to the hole along this 
accretion line. However, when pressure is included in the model, a cylindrical shock 
forms around this line and the accretion proceeds along an accretion column of 
high density and pressure shocked material. The predicted final accretion pattern 
consists of a stationary conical shock with the material inside the accretion radius 
being captured by the central object. An schematic representation of this solution 
is depicted in Fig. ^ 

A numerical evolution of relativistic wind accretion past a rapidly-rotating 
Kerr black hole (a ~ 0.999M, a specific angular momentum, M black hole mass) 
is depicted in Fig. (left panel). This simulation shows the steady-state pattern in 
the equatorial plane of the black hole. The tail shock appears stable to tangential 
oscillations, in contrast to Newtonian simulations with tiny accretors (see, e.g., B] 
and references there in; see Iq] for a related discussion). The accretion rates of mass 
and linear and angular momentum also show a stationary behavior (see, e.g., P, 0). 
As opposed to the non-rotating black hole, in the rotating case the shock becomes 
wrapped around the central accretor, the effect being more pronounced as the black 
hole angular momentum a increases. The inner boundary of the domain is located 
at r = l.Oilf {inside the event horizon which, for this model, is at 1.04A/) which is 
only possible with the adopted regular coordinate system. The fiow morphology 
shows smooth behavior when crossing the horizon, all matter fields being regular 
there. 
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Figure 2: Relativistic wind accretion on to a rapidly rotating Kerr black hole (a = 
0.999M, the hole spin is counter-clock wise) in Kerr-Schild coordinates (left panel). Initial 
model parameters: Uoo ~ 0.5, c^^ — 0.1 and F — 5/3. Isocontours of the logarithm of 
the density are plotted at the final stationary time t — 500M. The right panel shows 
how the flow solution looks like when transformed to Boyer-Lindquist coordinates. The 
shock appears here totally wrapped around the horizon of the black hole. The box is 
12M units long. The simulation employed a (r, (j!))-grid of 200 x 160 zones. 



The enhancement of the pressure in the post-shock zone is responsible for the 
"drag" force experienced by the accretor. The rotating black hole redistributes the 
high pressure area, with non-trivial effects on the nature of the drag force. The 
pressure enhancement is predominantly on the counter-rotating side. We observe 
a pressure difference of almost two orders of magnitude, along the axis normal 
to the asymptotic flow direction. The implication of this asymmetry is that a 
rotating hole moving across the interstellar medium (or accreting from a wind), 
will experience, on top of the drag force, a "lift" force, normal to its direction of 
motion (to the wind direction). Although different in origin this feature bears a 
superficial resemblance with the Magnus effect of classical fluid dynamics. 

The right panel of Fig. shows how the accretion pattern would look like 
were the computation performed using the more common (though singular) Boyer- 
Lindquist coordinates. The transformation induces a noticeable wrapping of the 
shock around the central hole. The shock would wrap infinitely many times before 
reaching the horizon. As a result, the computation in these coordinates would 
be much more challenging than in Kerr-Schild coordinates, particularly near the 
horizon. Since the last stable orbit approaches closely the horizon in the case of 
maximal rotation, the interesting scenario of co-rotating extreme Kerr accretion 
would be severely affected by the strong gradients which develop in the strong- 
field region. This will most certainly affect the accuracy and, potentially, also the 
stability of numerical codes. 
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